Transient simulation of regression rate on thrust regulation process in hybrid rocket motor

The main goal of this paper is to study the characteristics of regression rate of solid grain during thrust regulation process. For this purpose, an unsteady numerical model of regression rate is established. Gas–solid coupling is considered between the solid grain surface and combustion gas.Dynamic mesh is used to simulate the regression process of the solid fuel surface. Based on this model, numerical simulations on a H2O2/HTPB(hydroxyl-terminated polybutadiene) hybrid motor have been performed in the flow control process. The simulation results show that under the step change of the oxidizer mass flow rate condition, the regression rate cannot reach a stable value instantly because the flow field requires a short time period to adjust. The regression rate increases with the linear gain of oxidizer mass flow rate, and has a higher slope than the relative inlet function of oxidizer flow rate. A shorter regulation time can cause a higher regression rate during regulation process. The results also show that transient calculation can better simulate the instantaneous regression rate in the operation process.     

三维非对称斜劈高速流动计算方法研究

交叉激波是高超声速进气道发展的关键气动问题之一。为实现对此类高速流动的准确数值模拟,基于RANS方程的自研多块结构网格计算软件,对经典三维非对称斜劈构型的斜激波交叉干扰流动特征开展了计算研究。通过对比分析,系统评估了计算网格、湍流模型及高精度格式对空间流场、底板摩擦力流线、压强分布、绝热壁温分布等的影响,获得了适用于三维非对称高速流动问题的网格生成准则和计算策略,同时验证了自研软件的模拟能力。     

碳纤维复合材料基体改性/铜网组合雷击防护性能的研究

传统碳纤维复合材料（CFRP）树脂基体导电性差易遭受雷击损伤,本文使用石墨烯-镀镍碳纤维粉作为导电填料,对树脂基体进行电导率改性,并在表面铺设铜网,进行模拟雷电流冲击试验,检验基体改性/铜网组合雷击防护效果。试验结果表明,树脂基体改性后CFRP层压板在0°、90°纤维方向及厚度方向电导率分别为1.1571×10⁴、1.0871×10⁴、204.2 S/m,分别提高1.54倍、1.16倍、433.47倍。200 kA模拟雷电流A波冲击下,无防护试件雷击附着后明火燃烧,次生效应持续,而单一铜网防护和组合防护则能抑制次生效应;无防护表面最大损伤直径14.62 cm,此能量下铜网被击穿,单一铜网防护表面最大损伤直径19.05 cm,而组合防护表面最大损伤直径8.93 cm,下降53.12%;相比无防护试件,单一铜网和组合防护内部损伤面积分别下降66.2%和96.7%。单一铜网击穿后,树脂烧蚀后产生汽化反冲,增大损伤铜网脱落面积;组合防护铜网击穿后,改性树脂迅速导走电流,减小铜网脱落和内部烧蚀面积。     

基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

垂直/短矩起降飞机机翼内埋式风扇布局气动特性分析

以垂直/短距起降飞机过渡飞行状态为背景,针对机翼内埋式风扇布局的自由来流/风扇喷流混合型流动,基于结构/非结构混合网格使用CFD方法进行了非定常数值模拟和分析.首先使用滑移网格技术对NASA涵道螺旋桨进行算例验证,其时均计算结果与实验值的误差为5.3%,证明了计算方法的可靠性和准确性,然后数值模拟了机翼内埋式风扇布局在不同迎角下的气动性能.结果表明:风扇喷流在机翼上产生了特有的“抽吸”和“堵塞”效应,引起了机翼总升阻力的显著增加,升力最大增量达到干净机翼升力的2.6倍,阻力最大增量为干净机翼阻力的3.2倍,混合流场在机翼后缘引起了升力损失并卷起对涡.      

螺旋桨飞机滑流非定常数值模拟研究

采用非定常数值模拟方法,基于滑移网格技术,分析研究螺旋桨滑流对全机气动特性的干扰影响。用多块结构网格数值模拟软件,从搭接边界建立、搭接网格生成、流场插值传递3个方面对螺旋桨滑移搭接网格的数值模拟方法进行全面详细的介绍,并且阐述基于该技术的非定常数值模拟方法。在此基础上,开展螺旋桨飞机模型验证计算,分析研究螺旋桨滑流对全机流场和气动性能的影响,通过与风洞实验数据的对比分析,验证滑移搭接网格技术的可行性与准确性。研究结果表明,基于滑移搭接网格技术的螺旋桨滑流数值模拟,可以真实反映螺旋桨滑流流场特性,为螺旋桨类飞机的设计评估提供一种有效的数值模拟手段,在工程方面具有重要的实用价值。     

基于阶跃响应的飞行器单独静、动导数数值求解方法

基于阶跃响应方法,结合刚性动网格技术,对飞行器的单独静、动导数的精细化数值计算进行了研究.以纵向为例,通过给物面施加恒定附加攻角,求解得到阶跃响应运动过程的非定常气动力,求导得到静导数.同样给物面施加恒定的俯仰角速度,并同时强迫物面平动以抵消俯仰转动产生的附加攻角影响,可由非定常气动力求导得到动导数值.分别利用NACA0012翼型和三维SACCON飞翼无人机进行了计算验证,各攻角下的静、动导数值与文献、试验结果吻合得很好,最大误差不超过5%.结论表明:基于阶跃响应的单独静、动导数直接模拟方法计算耗时仅为传统强迫振动方法的21%,效率相对较高,且可推广到横航向的动导数计算,为飞行器的稳定性研究提供参考.      

基于在轨热环境的可展开天线反射面精度调整技术

随着宇航事业的发展,网状展开天线的尺寸越来越大,对形面精度要求越来越高(特别是高频段),在轨热环境与地面调整点的热环境差别大,必须采用基于在轨热环境的网面精度调整技术以满足反射面的精度要求。就在轨热环境下的辅助牵引面式网面精度调整技术的理论分析及调整方法进行了探讨,给出了计及在轨热环境的网面形状精度调整的基本理论公式及调整步骤,并对某样机进行了调整计算。计算结果表明,文中的分析、公式、方法是正确的,可行的。     

三角面元灰度图象三维重建算法研究

三角面元SFS算法是结合三角面元模型和线性反射图来阐述SFS问题的,该算法的基本思想是用三角面元素来近似光滑表面,并用一系列基于节点的线性方程组合来表达该表面。本文在此算法的基础上,求解时通过加权平均迭代,提高了恢复的准确性和迭代的收敛速度。     

高精度差分格式及多尺度流场特性的数值模拟

简要介绍了作者近年发展的高精度差分格式,给出了两类格式精度的分析方法,分析了粘性面差分逼近所引入的耗散量低亏效应。文中按群速度将格式分为快型,慢型和混合型三类。对每类格式分析了高频波的传播特性。